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ABSTRACT 

Kowalski & Saumon (2006) identified the missing absorption mechanism in the observed 
spectra of cool white dwarf stars as the Lya red wing formed by the collisions between 
atomic and molecular hydrogen and successfully explained entire spectra of many cool DA- 
type white dwarfs. Owing to the important astrophysical implications of this issue, we present 
here an independent assessment of the process. For this purpose, we compute free-free quasi- 
molecular absorption in Lyman-ff due to collisions with H and H2 within the one-perturber, 
quasi-static approximation. Line cross-sections are obtained using theoretical molecular po- 
tentials to describe the interaction between the radiating atom and the perturber. The variation 
of the electric-dipole transition moment with the interparticle distance is also considered. Six 
and two allowed electric dipole transitions due to H-H and H-H2 collisions, respectively, are 
taken into account. The new theoretical Lyman-ff line profiles are then incorporated in our 
stellar atmosphere program for the computation of synthetic spectra and colours of DA-type 
white dwarfs. Illustrative model atmospheres and spectral energy distributions are computed, 
which show that Lya broadening by atoms and molecules has a significant effect on the white 
dwarf atmosphere models. The inclusion of this collision-induced opacity significantly red- 
dens spectral energy distributions and affects the broadband colour indices for model atmo- 
spheres with T e g < 5000 K. These results confirm those previously obtained by Kowalski & 
Saumon (2006). Our study points out the need for reliable evaluations of H3 potential energy 
surfaces covering a large region of nuclear configurations, in order to obtain a better descrip- 
tion of H-H2 collisions and a more accurate evaluation of their influence on the spectrum of 
cool white dwarfs. 
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1 INTRODUCTION 

White dwarf (WD) stars represent the most common final stage 
of the stellar evolution, and as such convey valuable information 
about the history of our Galaxy and stellar populations (see Althaus 
et al. 2010 for a recent review). The large majority of observed 
WDs show hydrogen-rich atmospheres (DA WDs). Since the spec- 
trum of light escaping from old and cool WDs controls the rate at 
which they cool, the evaluation of cooling timescales of such WDs 
is sensitive to the gas/fluid opacity and equations of state used in 
the model atmospheres (Hansen 1998, Salaris et al. 2000, Serenelli 
et al. 2001). Furthermore, model atmospheres are a basic tool in 
the analysis of observed spectroscopy and photometric data. In this 
sense, the computation of synthetic spectra help to reveal important 
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details about the physical processes in external layer of WDs and 
to understand the spectral evolution of these objects. 

Photometric ob servations presented and analyzed by 
iBergeron et al.l i 19971) for 1 10 cool WDs showed a flux deficiency 
in the B magnitude for stars cooler than ss 5500 K. There, the 
UV absorption was interpreted to be a missing opacity source due 
to hydrogen. By computing the Lya quasi-molecular absorption 
and successful fi ts to the entire spe c tra of cool hydrogen atmo- 
spheres WDs, iKowalski & Saumo n (2006) identified the origin 
of B deficiency as the Lyman-a line broadened by collision of 
absorbing H atom with IF molecules. 

Quasi-molecular radiative transitions occur when a photon is 
absorbed or emitted by a hydrogen atom while it interacts with 
one or more neighboring particles (atom, ion or molecule) at 
atom-perturber separations of few Angstroms. The interaction be- 
tween absorber atom and other particle at short distances even- 



in the line wings 
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1982j). Lyman-a satellites have been observed 
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in the spectra of DA white dwarfs and identified as produced by 
H+H + (1405 A satellite ) and H+H (1623 A satell ite) collisions by 
iNelan & Wegneri Jl985h and iKoester et al.1 Jl985h . respectively. A 
broad H 2 collision induced satelli te in the red wing o f Lyman-/? at 
1 150 A has been also detected bv lAllard etal] d2004l) in the bright 
pulsating DA white dwarf G226-29. 

Classical resonance and van-der-Waals broadening have been 
demonstrated to be inadequate dSandoetal.lll969h in reproducing 
the Lyman a red wing extending far into the optical region. These 
broadening evaluations assume that the interaction between radiat- 
ing atom and perturber follows a power-law dependence r~ p with 
the interparticle distance r (e.g., p = 3 for resonance broadening 
and p = 6 for van-der-Waals broadening, c.f. Mihalas 1978), which 
fails to describe true interatomic potentials at r smaller than few 
Angstroms where far wing contributions are produced. Appropriate 
wing absorption theories are based on quasi-molecular approaches 
which use accurate theoretical molecular potentials to describe the 
interaction between the radiator and perturber and take into account 
the variation of electric-dipole transition moment with the inter- 
particle distance le.g.. Jablonski 1945; Bates 1951; Chen & Takeol 



Il957l ; iGallager & Holstemlll977l ; lAllard & Kielkopj|l982h . More 
specifically, the spectrum of absorption is evaluated using ab initio 
calculations of Born-Oppenheimer energies and dipole transition 
moments for the electronic states of the quasi-molecule or dimer 
formed during the collision. 

Atmosphere model evaluations for cool WDs that included 
Lyman a line br oadening based on qua si -molecule methods were 
performed by IKoester & Wolfti bOOd) . IWolff etalj J2002h and 
iKowalski & Saumorj J2006l) (hereafter KS). KS were successful in 
reproducing the spectr um of cool DA- WDs du e to the inclusion of 
broadening by H-H 2 . Ko ester & Wolff (2000) used only broaden- 
ing by H-H and H-He collisions. KS evaluations were based on the 
semi-classical approximation which is generally considered valid 
to describe the far wings of the line profile. This approach takes into 
accou nt the variation of the d ipole moment during an atomic colli- 
sion. [foesterj&Wolfll (2000)'s calculations have been made using 
the q uasi-static limit of the so-called unified theory l lAllard et al.l 
1 19991) . In the formalism of the unified theory, the whole line profile 
from the line core to the far wings can be assessed. It also takes 
into account additional contributions from multiple perturber colli- 
sion<[], which are considered important at perturber densities larger 
than 10 21 particles per cm 3 . The so-called static limit of the unified 
description produces a line spectrum similar to that of the semi- 
classical approximation. 

Given that more than 80% of the white dwarf population is 
of DA type (Eisenstein et al. 2006, Althaus et al. 2010), Lya wing 
broadening could have implications on the luminosity function of 
very old WDs. It is thus of much interest to study the wings of 
the Lyman-a absorption line yielded in atomic collisions. Consid- 
ering KS opacity calculations showed that the Lya quasi-molecular 
opacity may have an important role in affecting the emergent UV 
and blue radiation of cool DA- WD stars, we have performed an in- 
dependent evaluation. KS showed only results based on their es- 
timation of the main contributions to the gas opacity from H-H 
and H-H2 encounters. Here, we present results including all dipole- 
allowed transitions resulting from these collisions and show an 
overall analysis of the relative importance among them. For com- 
puting Lyo- wing opacity, we have elected the semi-classical broad- 
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ening theory as a matter of convenience at the present stage of our 
calculations, to reproduce KS calculation conditions and for com- 
parison purposes in future works. We postpone to develop wing 
broadening evaluations based on the unified method to a later work. 
As this paper is devoted to cool DA-type white dwarfs, we only in- 
clude collisions of radiating H atom with atoms and molecules. 

This paper is structured as follows. In Sect. [2] we briefly re- 
view the semi-classical broadening method and assess at some 
length the wing broadening arising from collisions of H atoms and 
H 2 molecule. The line opacity including non-ideal gas effects on 
the upper state of transitions are shown in Sect. 12.41 Model atmo- 
spheres and input physics are detailed in Sect. [3] We then evaluate 
in Sect. [4] the implications of our line broadening opacity on the 
emergent spectrum of cool white dwarfs. Finally, in Sect. [5] we 
summarize the main conclusions. 



2 LINE BROADENING THEORY 

In the electric dipole approximation, the cross section tr v of the 
radiation absorption corresponding to a transition from an initial 
state i to a final state j of an atom or molecule is given by 



o-ij(v) = 47r 2 av|r,/0, 7 (v), 



(1) 



where a is the fine-structure constant, v the frequency of the tran- 
sition, Tjj the electric-dipole transition moment, and 0y(v) the nor- 
malized line profile 



r 



(pij(v)dv= 1. 



(2) 



For radiation unpolarized and randomly oriented particles |r (J | 2 = 
\Djj(r), with Djj(r) the so-called dipole-strength function for the 
transition i — > /, As a reference, th e Lyman-a transition (i = Is, 
j = 2p) of isolated hydrogen atoms l lBethe & SalpeteJ 19571) . has a 
central frequency vo = 2.467 x 10 15 Hz, D\ Sy 2 P = (32/27) a 2 , (uq is 
the Bohr radius) and the resulting total cross section is 



f°° 4n 2 a 

0"is,2p(v)dv = ——v a D ls2p 

Jo J 



1.1044 x lQ- 2 cm 2 . 



(3) 



2.1 The quasi-static approach 

Our calculation of the wing absorption of Lyo- is made within the 
semi-classical approximation following the methodology used in 
KS. This approach is based on a number of assumptions 

(i) Born-Oppenheimer approximation: nuclei move on single 
adiabatic potential energy surfaces created by the much faster mov- 
ing electrons. 

(ii) Adiabatic collisions: the interaction between two particles is 
viewed as the formation of a quasi-molecule which moves over a 
particular Born-Oppenheimer energy curve during the collision. 

(iii) Nearest neighbor approximation: the radiating atom is as- 
sumed to be only perturbed by its nearest neighbor. 

(iv) (Classical) Franck-Condon principle: the radiative transi- 
tion occurs in the neighborhood of the internuclear distance where 
the difference between upper and lower potentials of the quasi- 
molecule equals the photon energy. 

The first enumerated approximation is used to identify collision- 
ally perturbed atomic states as parts of the initial and final molec- 
ular adiabatic states. The second assumption is considered valid 
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in slow atomic collisions ( Hilschfelde retaill 19541) and seems ap- 
propriate for the temperatures of cool WD atmospheres. In the 
one-perturber approach we ignore multiple simultaneous encoun- 
ters. Classically, the fourth assumption considers that the radiative 
decay time is short compared to the collision time. 

Within the preceeding assumptions, the quasi-static approx- 
imation predicts the following expression of the line profile 
dMargeneau & Lewis! 19591 ; I Afford & Kielkopj|l982h 

\dVij/dr\ 

Here, h is the Planck constant, and Vy(r) the difference between 
the two Bom-Oppenheimer energies Vj(r) and V,(r) representing 
the interaction of the active atom, in each one of two states i and j 
involved in the transition, with a perturber at a distance r. The first 
factor of the right hand side in Eq. {4]l takes into account the rela- 
tionship between atom-perturber distance and transition frequency 
y, which is given by 

hv = V u (r) = Vj(r) - V,(r). (5) 

The function Pi(r) in Eq. © is the probability density of finding 
the nearest perturber to a distance r from the radiator, 

P,(r) = Anrip? exp(-4mi p r 3 /3), (6) 

being n p the mean density of perturbers in the gas. The exponential 
factor in l[6]l is usually negligible numerically and is often omitted, 
however, it is necessary to obtain the correct normalization of the 
profile (pijiy) [Eq. Q]. 

With Eq. {4]l, the cross-section in Eq. (T) may be written as 



Table 1. H2-states asymptotically correlated with H(l.s)+H(ls, 2s, 2p). 



cr,/(v) 



An 1 a hv 
3 \dVi;/dr\ 



(7) 



where yS = (kT) , k is the Boltzmann constant and T the gas tem- 
perature. The exponential factor in Q is introduced as part of the 
Boltzmann distribution function of H atoms and accounts for the 
probability of finding perturbed atoms in the lower state i respect 
to unperturbed atoms [Vj(r — > oo) = 0]. If the potential difference 
Vjj(r) has an irregular behavior and there are several distances r ver- 
ifying the relationship {5), we have to sum over all atom-perturber 
configurations that contribute to the wing-broadening at the fre- 
quency v. 

The total cross section of a specific collision-induced absorp- 
tion line is a sum over all lower (i) and upper (J) molecular states 
which can contribute to the wing broadening. These contributions 
are additive in the Born-Oppenheimer approximation with appro- 
priate statistical weights Try. In particular, there are many contri- 
butions to the resonance broadening of the Lyman-or line corre- 
sponding to transitions of different quasi-molecular aggregates. In 
cool DA white dwarfs, the main perturbers are H atoms and IT 
molecules because they are the most abundant species. In the fol- 
lowing sections we consider wing broadening arising from colli- 
sions of H atoms with these particles. 

2.2 H-H collisions 

Table [T] lists the products of H(Ls)-H(ls, 2s, 2p) adiabatic colli- 
sions, including the equilibrium electronic energy T e for ground 

2 This approximation breaks down for high-energy collisions which can 
result in either ionization or excitation to other electronic states. Such col- 
lisions constitute non-adiabatic processes and are considered negligible in 
the present work. 
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rotovibrational levels as given bv lField etal]jl966h . the multiplic- 
ity g e of each state, and the data source for energy curves relevant 
to the present work. 

According to the adiabatic approximation, from the encounter 
of two H atoms, one in the state 1 s and the other in the state 2p, one 
of eight possible H 2 electronic states is formed (Table[TJ. However, 
only four (S'E„, a 3 E g , Cn„ and i 3 Tl g ) of these molecular levels 
have allowed electric dipole transitions with the lowest H2 states, 
X l ?.g and b 3 !,*, which asymptotically correlate with two sepa- 
rated H(li') atoms. Similarly, there are four H 2 levels which can 
be formed from adiabatic collisions of H(l j) and U(2s) atoms (Ta- 
bleQJ, but only two (fi''E„ and h 3 T. g ) of these states have allowed 
dipolar transitions to states X'Sg and £> 3 X^. Only states involved in 
molecular transitions allowed by electronic selection rules are of 
interest in the present study. 

The complete set of electric-dipole radiative processes con- 
tributing to hya wing from H-H collisions is illustrated in Fig. [T] 
The conditional probability P(i\i — > j), with values shown in Fig.[T] 
is the probability that a quasi-molecule in the state i is formed from 
an specific H(l.s)-H(Ls, 2s or 2p) collision, under the condition that 
i — » j is an allowed electric-dipole transition contributing to the 
Lya- wing (similar definition follows for P(j\j — > 0)- The classical 
evaluation of P(i\i — > j) takes into account that if two particles ap- 
proach adiabatically each other, the fraction of occasions on which 
they move along a particular energy curve (says z) is given by the 
ratio of the statistical weight of this curve (g/) to the sum of the 
statistical weights of all possible curves. 

The Born-Oppenheimer energy curves arising from H-H col- 
lisions that contribute to the Lya- wings are shown in Fig. [2] The 
potential curves adopted in this work are referenced in TableQ] As 
the atom and perturber get further away from each other, the elec- 
tronic energies tend to asymptotic values, which are sums of in- 
dividual particle energies. Since atomic states are degenerate, there 
are in general several molecular energy curves that tend to the same 
asymptotic energy. 

The dipole moments of molecular transitions allowed by elec- 
tronic selection rules are shown in Fig. [3] The dipole-forbidden 
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collision i-state j-state collision 

P(ili-j) P(jlj-i) 



H ls +H ls 




.b 3 E„ 



X'2. 




l/g^ H ls + H 2p 



h 3 E g >s^V 4 



H ls +H 2s 



Figure 1. Collision diagram for H(ls)+H(ls, 2s, 2p) encounters with inter- 
mediate quasi-molecular states connected by allowed electric-dipole transi- 
tions (dashed lines). Conditional probabilities of electronic state formation 
(see text) are indicated on the plot. 



10 



> 

CD 

>> 
OD 
U 
CD 

a 
w 





n 1 i 1 1 1 i 

BA i3n g 


1 1 1 1 1 1 

H(ls) + H(n = 2) - 
- 












B»S U 






H(ls) + H(ls) " 
>■ 




X'E 

1 ' 
* / 
\ * 






\ * 

\S 

1 1 1 , 1 1 1 





2 4 6 8 

r [bohr] 

Figure 2. Born-Oppenheimer energies of H 2 states relevant to this paper. 
Dashed lines are for singlet X states, solid lines for triplet £ states, long 
dashed line for the C IIJ state and the dotted line for the i 3 Yl* state. The 
asymptotic states corresponding to great internuclear separations are also 
indicated. 



character of the transition (Is - 2s) for separated atoms, is associ- 
ated to electric-dipole transition moments of the H 2 molecule T X -b- 
and T b _ h which become zero as the internuclear separations r — > 00. 
On the other hand, the correct dipole transition moment for the 
\s - 2p atomic transition is obtained in the dissociation limit for 
X — B,X - C, a - b and h - b transitions. 
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Figure 3. |7j;| 2 dipole moments between H 2 -states which correlates with 
H(2s or 2p) + HQ s) and H(ls) + H(ls) as a function of the internuclear 
distance. Data from|PressleiMl985t) for X -Ban d X - C. lFord et all Jl975l) 
for X-B'. and lStaszewska & Wolniewiczl ( 1999) for triplet state transitions. 



To obtain the total absorption coefficient cr H - H (v) due to H-H 
encounters, a sum of ifTJ over all H2 transitions which contribute at 
the frequency v must be performed, 



(Th-h(v) = ^ijCTijiv). 



(8) 



In the present case, the weights n\j of these contributions are di- 
rectly related to the probabilities P(j\j — * ')• It can be demonstrated 
that© recovers the correct limit value for asymptotically separated 
atoms as given by {3). 



2.3 H-H 2 collisions 

Adiabatic collisions between H atoms and H 2 molecules form tri- 
atomic systems. The analysis for atom-atom collisions in Section 
I2.2l was simplified by the fact that only a single internuclear coordi- 
nate, the atom-atom distance r, need be considered. In the present 
case, however, collision processes involving H 2 molecule and H 
atoms are functions of the interparticle distance r and also of the 
molecular orientation, which can be characterized by a single an- 
gle 9 in the case of homopolar diatomic molecule (9 is the angle 
defined between the line connecting the H atom to the bisector of 
molecular axis, see geometrical configuration insert in Fig.|5j- Con- 
sequently, the Bohr-Oppenheimer energy solutions for the poly- 
atomic molecule H 3 consist on energy surfaces instead of energy 

curves. 

Following IPeng et alJ l l 19951) . we denote E\, £ 2 , £3 and £4 
the four lowest potential energy surfaces of the H3 system. E[ la- 
bels the ground electronic state which asymptotically separates into 
H 2 (X'X g )+H(li). The first excited state of H3, E%, adiabatically dis- 
sociates into H 2 (fo 3 E„)+H(l.s) and undergoes a conical intersection 
with the ground state at equilateral triangle (D ih ) geometries. E\ 
and E2 become the degenerate 2p 2 E' state at Dy, structures and 
constitute a well-known Jahn-Teller system. The next two excited 
states, E 3 and E4 (2s 2 A\ and 2p 2 A'^ states in D 3h symmetry), cor- 
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Figure 4. Electronic states of H3 which in the asymptotic region H+H2 cor- 
relate with H(2s or 2p) + H 2 (X 1 J.+) and H(ls) + H 2 (X'X+) as a function of 
the H-H2 inte rparticle distance r and for co llision angle s 8 = 45° and 90°. 
Data f rom KS ilKowalski & Saumorj j200d) (fig. 1 ) KG:lKulander & Guesil 
\ 19791) (fig.D. lRoach & Kuntzlll986l) (fig. 5) Ma:lMavne et al.Nl984l) (an- 
alytic fit), and TH: lTruhlar & Horowitzl U978, 1979) (analytic fit). 



relates with H 2 (X%)+H(2s) and H 2 (X%)+H(2p), respectively. 
These states are the lowest Rydberg states of H3 and are located 
very close in energy for all geometries. 

It must be noted that the (£2) low-lying excited state is not 
relevant to the present study because it correlates with the b 3 l.„ un- 
bound electronic state of the H2-molecule rather than excited states 
of the H-atorr0. Therefore, the polyatomic states of interest here 
are £1, £3 and £4. These H3 states are sufficient to infer the likely 
course of an encounter between H and H2 contributing to the Lya 
wing absorption (Mayne et al. 1984). 

The evaluation of potential energy surfaces (PES) is com- 
putationally expensive compared with similar studies of diatomic 
molecules. Theoretical calculations are often depending on the 
quantum chemical method applied and, in general, it is not easy 
to identify the best available PES, specially for excited electronic 
states. Moreover, in the case of triatomic hydrogen, excited state 
energies are often available for a few nuclear configurations so 
that interpolations or fit procedures are required to cover a large 
region of the nuclear geometry as it is necessary in t he analysis 
of H-H 2 collisions. Accurate e nergy data obtained by iLiul d 19731) 
and by Siegbah n & Lid (fl978l) for the H3 ground state have been 
fitted by iTruhlar & Horowit j dl978l Il979h . This fit is considered 
one of the four most accu rate and widely used PES for the H 3 
system dMielke et alj |2002). Unfortunately, ener gy data of the first 
Rydbe rg states are limited to a few geometries. iKulander & Guesj 
l l 19791) car ried out calculat i ons al ong an equilateral insertion path 
(0 = 90°). lRoach & Kuntzl dl986l) used the diatomic-in-molecules 
(DIM) procedure to calculate the ground and lower-lying excited 
PESs over several geometries, p articularly for a 45° a ngle of ap- 
proach. Energy computations of Pets alakis et al.l dl988l) were per- 



The Ej state could be involved in ternary collisions since it asymptoti- 
cally correlates with 3H(ls). 



formed mainly for collinear H+H 2 approach (0 = 0°) and linear 
symmetric H3, and for a few perpendicul ar trajectories. Mo re recent 
evaluations of Rydberg energies due to IPeng et al] d 19951) contain 
also few data points. A s emiempirical PES o f the first H3 Rydberg 
states was developed bv lMavne et alj 1 1984) based on e valuations 
in the DIM approximation (Ravnor & Herschbach 1982). 

Some of these energy evaluations are shown in Fig.|4]for non 
collinear configurations, = 90° being the more probable condition 
of impact. Truhlar & Horowitz (TH) analytical representation for 
the ground state agrees remarkably well with DIM evaluations of 
Roach & Kuntz (RK) at = 45°, and with theoretical calculations 
of Kulander & Guest (KG) at = 90°. The equilibrium internuclear 
separation of the H2 molecule was fixed at R = 1.4 Bohr (0.74 
A), except in the case of KG evaluations which are based on the 
equilibrium distance of the H3 molecule (R = 1.65 Bohr, this is 
only slightly larger than the bond length of 1. 62 Bohr for H 3 ). Fig . 
[4] also shows values of the analytic fit due to iMavne et al.l d 19841) 
(Ma) for the first Rydberg states at = 90° (equivalent energy is 
assumed for £3 and £4 in this work). Ma fit yields smaller values 
than those from KG calculations throughout the region 2 < r < 4.5. 
At H-H 2 separations w 2 Bohr, the Ma curve crosses the KG curve 
and remains above of that for shorter distances. 

For comparison we have a lso plotted in Fig.|4]the results used 
b vlKowalski & Saumonl (2006) (KS), based on theoretical studies 
of lBoothrovd et al.l a 199 lL 1 19960 . We have not incorporated results 
of Boothroyd et al. in our wing opacity evaluations because no in- 
formation is available to which asymptotic state the evaluated lev- 
els correlate with, and there are not enough analyzed geometrical 
configurations to correctly identify the £3 and £ 4 states from the 
computed data points (Boothroyd 2008, private communication). 
The ground energy curve used by KS at = 90° lies above the TH 
and KG evaluations at H-H 2 separations smaller than 1.7 Bohr. For 
the first Rydberg state, KS results remain above those of KG at dis- 
tances smaller than about 5 Bohr with deviations of =s 0.8 eV for 
intermediate distances up to 1.5 Bohr and increasing differences at 
smaller distances. 

From this limited amount of PESs, one is already able to see 
the complexity in the selection of appropriate PESs to be used in 
opacity evaluations. Three important facts deserve to be mentioned: 
(0 r < 2 Bohr represents the critical region where absorptions to 
optical spectrum could be yielded, (») energy differences between 
E[ and £ 2 , £3 states play a decisive role in defining this spectrum, 
and (m), among the available PES data, KG and RK results provide 
simultaneous and self-consistent evaluations of energies for these 
states. Hence, we decided to adopt KG and RK data to develop 
PESs for the ground and lowest Rydberg states. In practice, we have 
considered the £3 and £4 potential energies identical because they 
are similar for all molecular geometries (differences smaller than 
1000 cirr 1 , May ne et al.ldl984l) ) and because they cannot be clearly 
distinguished from one another in the available data. In order to 
construct PESs covering a large range of nuclear geometries, we 
use angular interpolation between KG (0 = 90°) and RK (0 = 45°) 
evaluations. In this sense, it is worth noting that the largest con- 
tribution to the absorption cross section comes from perpendicular 
collision paths (the averaged-angle cross section, see Eq. (9} later in 
this Section, concentrates the integration weight near the configura- 
tion = 90°). To complete the PESs, we maintained the RK values 
for < 45°. We have verified that this choice introduces a minor 
uncertainty in the cross section computation. The consequences for 
the wing absorption calculation due to the adopted PESs and the 
use of alternative PESs based on TH and Ma fits will be discussed 
in Sect. |4] 



6 R. D. Rohrmann, L. G. Althaus and S. O. Kepler 



0.8 

J3 0.6 
o 

==. 0.4 
0.2 




- 




— i 1 t r r 1 1 r 1 

1 1 1 1 1 1 1 1 1 1 1 1 

6=0'(triangles) 






6=45*(crosses) 


_ 




0=9O*(circles) 




\ El 








H - 










/ * \ 


H« — - 

\ R 


1 1 1 


™ i i i 1 i i i 


1 i i i 1 i i i 1 i i i 



2 4 6 8 10 12 
r [bohr] 



priate for atmosphere model calculations. These results are avail- 
able in the web site http://www.fcaglp.unlp.edu.ar/evolgroup 

The total Lyor wing profile cr, or (v) is sensitive to the relative 
abundance of the different perturbers responsible for the line broad- 
ening. Within the nearest neighbor quasi-static approach, the rela- 
tive contributions to the line broadening by H and H 2 are given 
by the probabilities of finding an atom or molecule as the clos- 
est neighbor of radiating atoms. These probabilities can be approx- 
imated by the molar fractions of atoms (x H ) and molecule (xn 2 ) 
de.g. Rohrmann & Zoredl2006l) . resulting in 

o"m(v) = Xh(Th-h{v) + x H2 a- H -H 2 (y). (10) 

Finally, the extinction coefficient (units of cirT 1 ) due to collisional 
induced absorptions in the Lyman-o- line is 

X(v) = n H( „ =i) a m {v), (11) 

where n H ( n =i) is the number density (in cnT 3 ) of atoms in the 
ground state. 



Figure 5. |7y| dipole moments of Ei —E3 and E\ -£4 polyatomic transitions 
as a function of the nuclear geometry. Data from Petsalakis et al. (1988) and 
Peng et al. (1995) with fill and open symbols, respectively. Lines indicate 
our fits to results of Petsalakis et al. and Peng et al. for = 0°. 

Although the energy surfaces for the £3 and £ 4 excited elec- 
tronic states are nearly identical, the electric dipole moment be- 
tween these states and the ground state differs significantly. In the 
current s tudy, we adopt dipole m om ent transitions bas ed on evalu- 
ations of IPetsalakis et alj dl988h and lPeng etal] Jl995h . Analytical 
fits including internuclear distances and collision angles beyond the 
data range were used in our calculations. Fig.[5]shows that the \T t j\ 
transition moment between the E\ and £ 4 states increases with the 
atom-diatomic distance from r ~ 1.8 Bohr and approaches its the- 
oretical value of 0.74 Bohr when r — » 00. In D ih symmetry config- 
uration this transition is forbidden. At short interparticle distances 
(r < 4 Bohr), the transition from the ground state to the £3 state is 
predicted to be considerably stronger than the transition to the £ 4 
state. Although the £1 - £3 transition is allowed, it decreases above 
r ~ 2 Bohr and approaches zero for r — > 00, in agreement with 
the expected in the dissociation limit where the Is — > 2s atomic 
transition is forbidden. 

Finally, the total cross-section for Lyo- wing broadening by 
H-H 2 collisions is averaged over all collision-frames angles and 
expressed as 

o-h-h 2 (v) = X f o-ij(y,e)smed6 (9) 

where the summation comprises E t — > £3 and £j — » £ 4 transitions 
and a-, 7 (v, 0) is given by Eq. @ with V i} (r, 9), D y (r, ff) and V;(r, ff) 
functions of both the distance r and angle orientation 6 of the dimer. 
The Ttjj weights were chosen so that, using the asymptotic values of 
H 3 transition moment, the isolate atom limit is recovered [Eq. {3)]. 

2.4 Total collisional-induced wing absorption 

The lines profiles cr H - H (v) and a H -H 2 {v) given by Eqs. <[8j and l[9j 
have been convolved with a Doppler profile to take into account the 
thermal broadening introduced by the particle motions. We have 
prepared detailed opacity tabulations of these cross-sections appro- 



3 MODEL ATMOSPHERES 

The LTE model atmosphere code used in our analysi s is a mod- 
ified version of t hat de scribed at length in iRohrmannl d200lh and 
iRohrmann et al.1 J2002[) , which is appropriate for hydrogen and he- 
lium atmospheric compositions (including mixed and pure mod- 
els). Models are computed assuming hydrostatic and radiative- 
convective equilibrium. Energy transport by convection present in 
the cool atmospheres here considered is treated within the usual 
mixing length theory, where we have assumed the so-called ML2 
parameterization of the convective flux. 

The gas model used in the code includes the species H, IT, 
H + , IF, H+, H+, He, He , He + , He 2+ , HeJ, HeH and free electrons. 
The relative abundances of th ese species is determined b y the occu- 
pation probability formalism dHummer & M ihalas 1988). The cal- 
culated level occupation probabilities are then explicitly included 
in the calculation of the line and continuum opacities as described 
in Section l3TI 

The Lyo' wing opacity analyzed in the preceeding sections 
represents a partial contribution to the total gas opacity in atmo- 
spheres of cool DA WDs. For pure hydrogen models, the opacity 
sources in the numerical code include also bound-free (H, H~, HJ) 
and free-free (H, H 2 , H 3 , H~, H^, HJ) transitions, electronic and 
Rayleig h (H, H 2 ) sc attering, collision-induced absorptions (CIA) 
of H-H toovldfl968l). H-H, dGustafsson & Frommholdll2003l) and 
H 2 -H 2 dBorvsow et au 2001). and the most significant H line series. 
Details are described at long in lRohrmann et al] d2002l) . It should 
be pointed out that the H 2 -H 2 CIA consists on rotovibrational tran- 
sitions of H 2 induced by molecular collisions, which has a strong 
influence on the infrared spectrum of cool WDs. 

3.1 Nonideal effects in opacity laws 

One important advance in white dwarf atmosphere modeling has 
been provi ded by the use of the occu pation probability formal- 
ism due to Hummer & Mihalas (1988) (HM). The HM approach 
considers the perturbations on each atom or molecule by charged 
and neutral particles and includes their effects in the evaluation of 
atomic populations and equations of state of the gas. The interac- 
tions with neutral particles are treated within the hard sphere model 
and those from charged particles are calculated with micro-field 
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Figure 6. Occupation probabilities of the n = 1 (dotted line) and n = 2 
(dashed line) hydrogen states, and of the ground H2 molecule state (long 
dashed line) calculated in the photosphere (tr ojj = 1) of logg = 8, hydro- 
gen pure models as a function of effective temperature. 

distribution functions. The internal partition function of a given 
specie is then written as 



(12) 



where E t and g t are, respectively, the excitation energy and mul- 
tiplicity of the level i. The function w, is the so-called probability 
occupational of the level and it is computed self-consistently with 
the non-ideal term in the gas free energy. The w, values decreases 
continuously and monotonically as the strength of the relevant in- 
teraction increases, and avoid the familiar divergence of internal 
partition functions. Analytical continuity of all thermodynamical 
properties of the gas is assured by the application of the free energy 
minimization tec hnique. Details of our HM evaluations are given in 
iRohrmann et alj I2OO2I) . Fig. [6] shows the occupation probabilities 
of the lowest bound states of H and IT at a Rosseland mean optical 
depth t Ross = 1 for model atmospheres at 2400 K < T eff < 8000 K 
and logg = 8. Considering that the dominant perturbers in atmo- 
spheres of cool WDs are neutral particles, the occupation probabili- 
ties are mainly determined by the hard sphere model. The transition 
of Wi from near unity to near zero is a monotonically decreasing 
function of the gas density which, in Fig. [6] increases with decreas- 
ing r cff . 

Because the pressure s hifts observed exp erimentally for hy- 
drogen lines are very small l lWiese et al.ll 19721 ) and because of the 
lack of a reliable theory to compute energy level shifts in bound 
particle states, the HM formalism uses energy eigenvalues of iso- 
lated particle. Non-ideal effects are therefore directly accounted for 
in Wj factors leading to the concept of an effective statistical weight 
for bound states, whereas the internal particle structure is assumed 
unperturbed. This includes the use of oscillator strengths of isolated 
atoms. However, since the relative population between two particle 
levels is modified respect to ideal (Boltzmann or Saha) relations [as 
it can be inferred from Eq. J12H. the atomic transition rates must be 
accordingly adapted to satisfy the principle of detailed balancing at 
thermodynamic equilibrium. If we denote Ay as the transition prob- 
ability per second from ; to j levels of an unperturbed atom, then 



tlihijWj. 



(13) 



This phenomenological proposal preserves the well-known Ein- 
stein and Einstein-Milne relations on the Ay coefficients. Thus, 
AjjWj plays the role of a conditional probability such that it ap- 
proaches to Ay (the value corresponding to isolated particles) for 
non-perturbed transition final states, and approaches to zero for 
strongly perturbed levels. Expression d!3t is applied to all bound- 
bound and bound-free transitions with the convention that bare ions 
have formally Wj = 1 . 

The evaluation of Lyce wing opacity (and CIA processes 
in general) departs from the non-ideal opacity laws of the HM 
method. This happens because quasi-molecular absorptions arise 
from processes that involve atoms strongly disturbed. In fact, parti- 
cle interactions which include variations of both energy and dipole 
moment of the transition need be directl y considered in th e evalu- 
ation of such absorption cross sections JKowalski 2006b). There- 
fore, it would be incorrect to use the non-ideal opacity law H31 in 
the context of far-red wing absorption since particles perturbations 
are already explicitly included in expression dlOt - Of course, HM 
approach participates in the absorption coefficient il lb as it decides 
the actual populations of ground state atoms and perturbers which 
take part in the line wing broadening. 



4 RESULTS 

The properties of the hya broadening can be understood by study- 
ing Figs. UJ and [8] which show the absorption cross section of 
Lycr line broadened by collisions with H and H 2 , respectively, for 
T = 6000 K and perturber density n p = 10 2C cm , These figures 
illustrate also the photon wavelength for transitions induced in H- 
H and H-H 2 encounters according to the classical Franck-Condon 
principle [Eq. {5}]. 

For H-H collisions, the difference between the upper and 
lower interatomic potential for b — a,b — i and b — h triplet transi- 
tions increases (and the wavelength decreases) monotonically with 
the internuclear separation (bottom panel of Fig.UJ. Consequently, 
their wing spectrum contributions are relatively featureless and ex- 
tend to long wavelengths from the Lycr central wavelength at 121.6 
nm (top panel of Fig.UJ. The potential energy difference for X — B, 
X-C and X - B' transitions does not decrease monotonically with 
the internuclear separation. Local extrema in the X - B, X - C and 
X — B' energy differences correspond respectively to wavelengths 
A = 162.3, 126.8, and 121.8 nm at the red wing, and A = 83.0, 82.0 
and 73.4 nm at the blue wing. Each extrema can eventually pro- 
duce a satelli te feature in t he line wings (a satellite at 162.3 nm was 
pr edicted by Sando et al. ( 1969) and observed in a WD spectrum 
bv lKoesteretal.Ul985h . These energy extrema yield the classical 
discontinuities at wavelengths observed in the absorption cross sec- 
tions for singlet transitions (Fig. [7}- The blue wing due to singlet 
transitions is caused by the deep well in the ground electronic state 
of the diatomic molecule. In the far red wing, the b - a transition 
gives the dominant contribution to the line broadening by H-H col- 
lisions. Comparison with results of Kowalski & Saumon (symbols 



4 Notice that the proposal of transitions from bou nd states to the so-called 
dissolved states considered in Hubeny et al. 1 1994) is invalid within of HM 
formalism. See details in Ro hrmann et alj 1 2002). 
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Figure 7. Above: Cross-sections as given by the quasi-static approach, Eq. 
(7), for different transitions in H-H collisions. Dashed line for -fi'Xj, 
long dashed line forX 1 X+-C 1 n+, and dotted line for Results 
due to transitions between triplet states are labeled on the plots (dot-dashed 
lines). The density of perturbers (H atoms) is n p = 10 20 cnT 3 and the gas 
temperature T = 6000 K. Cross section evaluations of Kowalski (2006, fig. 
49) are shown with symbols. Below: Variation of the wavelength based on 
the level energy difference [Eq. {5}] of selected transitions for H-H dimers 
as a function of the internuclear distance. The curves for b - i and b — h 
transitions are visually identical. 



Figure 8. Above: Cross-sections as given by the quasi-static approach, Eq. 
{7), for E[ - E3 and E\ - £4 transitions in H 2 -H collisions. Notation in the 
text. The density of perturbers (H2 molecules) is n p = 10 20 cm and the 
gas temperature T = 6000 K. KS results are taken from Kowalski (2006, fig. 
49). Below: Variation of the wavelength based on the level energy difference 
[Eq. (3}] of select H3 transitions as a function of the internuclear distance 
and for different PES evaluations (indicated on the plot). KG-RK curves 
correspond to 9 = 45°, 60°, 75° and 90°. KS results are based on fig. 48 of 
Kowalski (2006). 



in top panel of Fig.O shows very good agreement for the dominant 
absor ption at the far red wing. 

Allard & Kielkopf (2009) noted in a recent study that simul- 
taneous collisions with more than one perturber gives an additional 
absorption from X-B transitions. Specifically, collisions with mul- 
tiple atoms yield X-B opacity contributions which dominates 
the 180 - 300 nm region. These multiple collisions could intro- 
duce additional absorption features in the UV spectra of DA WDs, 
however, they have likely no impact over the UBVRI photometry 
of these stars because the corresponding bandpasses are at longer 
wavelengths. Furthermore, multiple perturber effects become im- 
portant at densities as high as 10 21 particles per cm 3 , but the atomic 
population hardly reaches this value in WD atmospheres due to 
molecular recombination. 

Top panel of Fig. [8] shows the absorption coefficients from 
Ei — > £3 and Ei — > £ 4 transitions in H-H 2 collisions as a func- 
tion of the wavelength for T = 6000 K and molecular density n Hl = 
10 20 cnT 3 . The variations of the wavelength with the H-H2 distance 
is also given in Fig.[8](bottom panel). KG-RK la bel denotes results 
based on i nterpolations of PES calculated by iKulander & Guesq 
jl979h and lRoach & Kuntj jl986l) as described in Sect. [23] while 
Ma-TH label corresponds to results obtained from analytical fi ts of 
Ei and £ 3 (fi) PES s due to lTruhlar & Horowitz! Jl978lll979h and 
iMavne et alj dl984l) . respectively. 

According to the semiclassical approximation, no broadening 
of the Lyo- transition due to H 2 molecules will occur for A < 121.6 
nm (bottom panel in Fig. [H}. Cross sections based on KG-RK en- 
ergy surfaces (solid lines in top panel of Fig. [8} show that the far 



wing broadening by H9 is dominated by the transition Ei - £3 and 
originated from short range interactions (r < 2 Bohr). The E t - £4 
transition gives an important contribution in the near wing (A < 150 
nm) but yields smaller cross section values than £1 -£3 transition to 
longer wavelengths. As a remarkable feature, the cut-off at A = 400 
nm in the £1 - £4 absorption is caused by the reduction of the dipo- 
lar moment at short internuclear distances (Fig. [5}. 

Fig.[8](bottom panel) shows that the wavelengths derived from 
KS potentials at - 90° are below the KG-RK values at relatively 
great distances (1.5-4 Bohr), above it at intermediate distances (1.3- 
1.5 Bohr), and finally falls at short distances (r < 1.3 Bohr). Such 
differences are relatively small in H-H 2 separations where the far 
wing is formed so that both calculations yield similar cross sec- 
tions, as it can be appreciated in Fig.[8](top panel). Compared with 
results based on KG-RK energies, KS absorptions increase slightly 
towards long wavelengths (A > 450 nm) likely as a result of the 
energy discrepancies found between 1.3 and 1.5 Bohr. 

Fig. [8] also shows that the opacities calculated with Ma-TH 
PESs (dotted lines) are similar to other reported results for wave- 
lengths shorter than 250 nm but then fall sharply at longer wave- 
lengths. Discrepancies are particularly severe for the E\ - £3 tran- 
sition in the optical region (A > 400 nm). In this figure (bottom 
panel), one can see that the cross section at visible wavelengths is 
governed primarily by the energy difference between initial and fi- 
nal states in the H 3 transitions at atom-molecule separations smaller 
than 2 Bohr. At these distances, wavelengths derived from Ma- 
TH PESs noticeably differ from KG-RK and KS results. In par- 
ticular, Ma-TH curves exhibit a local maximum about 400 nm 
for 9 = 90° (shorter wavelengths for decreasing collision angles). 
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Figure 9. Monocromathic opacity coefficient in the photosphere (tr oss = 1 ) 
of H model atmospheres at log g = 8 and r c fj = 7000 K (upper), 5000 
K (medium) and 3000 K (lower). Shown are the total continuum opacity 
without hya absorption (dashed lines), the individual contributions to hya 
wing broadening by H (dot-dashed lin es) and H2 (solid li ne), and the Stark- 
broadened hya profile (dotted lines) o f lVidal et alii 19731) . Temperature and 
gas density values are indicated on the plot. 



Figure 10. Emergent spectra calculated for hydrogen pure atmosphere mod- 
els with (solid lines) and without (dashed line) collision broadening of hya 
at logg = 8 and 7W = 7000 K (upper), 5000 K (medium) and 3000 K 
(lower). Dotted lines represent black body spectra at T = T t g. For compar- 
ison, are shown with symbols the synthetic spectra computed by Kowalski 
(2007) for T eff = 3000 and 5000 K. Short solid lines indicate the location 
of the transmission functions for the filters F606W and F8 1 8W of the HST 
ACS (Vega-mag system) and those of the U BVR1 photometry. 



Consequently, opacities based on Ma-TH energies have a dramatic 
long-wavelength reduction as shown in the top panel of Fig. [8] Fur- 
ther analysis shows that these discrepancies are mainly due to the 
different behaviors of the PESs used for the Rydberg states. 

Taking into account that absorption cross section values are 
very sensitive to small changes in the difference between the PESs 
of H 3 states, we are inclined to believe that KG and RK evalua- 
tions (which give simultaneous evaluations of E t and E 3 energies) 
provide a better internal consistency than the use of TH and Ma an- 
alytical fits, which are based on different energy computations for 
Ei and £ 3 states. A conclusion of our study is that reliable evalua- 
tions of H 3 PESs covering a large region of nuclear configurations 
are needed for a better description of H-H2 collisions and an accu- 
rate evaluation of their influence on Lya wing broadening. Results 
shown henceforth are based on KG-RK energies. 

The importance of the collision induced broadening of hya 
line can be appreciated in Fig. [9] which displays the total continu- 
ous monochromatic opacity (dashed lines) and the Lya wing broad- 
ening by H (dot-dashed lines) and H 2 (solid lines) as a function 
of the wavelength for different physical conditions of a hydrogen 
gas. These evaluations correspond to the Rosseland optical depth 
T Ross = 1 of pure hydrogen model atmospheres at logg = 8 and 
r cff = 7000 K (upper), 5000 K (medium) and 3000 K (lower). For 
the sake of compa rison, we depict the Stark broadening Lya profile 
( Vidal et al.ll 19731) (dotted lines), that accounts for the main source 
of broadening in most relatively hot WDs (T c g > 10000 K). This 
broadening mechanism, which is due to charged particle interac- 
tions, markedly weakens for the low T cff analyzed here. 

For the highest temperature considered in Fig. [9] T = 7541 
K, the absorptions originated from singlet transitions in H-H in- 
teractions dominates the gas opacity coefficient between 91 nm 



(the Lyman jump) and 160 nm. The red wing of Lya is sensi- 
tive to the degree of molecular recombination in the gas, because 
that determines the relative importance of broadening by atoms and 
molecules. At temperatures lower than T < 5000 K and densities 
greater than p « 0.0003 g cm" 3 , the degree of dissociation has been 
reduced to 40% and the far red wing exhibits an enhanced contribu- 
tion due to H-H 2 collisions. Below T = 4000 K, the atoms are sur- 
rounded mostly by molecules and the broadening by H 2 dominates 
the total gas opacity from the Lya core up to about 500 nm, while 
H-H absorption remains as the main opacity source at wavelengths 
blueward of the Lya line. These results ba sed on KG-RK ener- 
gies ar e in agreement with those obtained by Kowalsk i & Saumoril 
d2006h . 

Synthetic spectra of DA- WD atmospheres with logg = 8 and 
r eff = 3000, 5000 and 7000 K are show in Fig.[l0] The short wave- 
length spectrum of the = 7000 K model is clearly affected 
by H 2 quasi-molecular absorption between Lya core and « 200 
nm. The only prominent feature in this spectral region is the lit- 
tle jump at 160nm originated by the - transition. At 
r cff = 5000 K, the contribution from H 3 transitions becomes domi- 
nant far away of the line center as described above and significantly 
weakens the emergent radiation for A < 450 nm. At the coolest 
model, r cff = 3000 K, the broadening by H 2 results in a Lya wing 
extending far into the optical region (up to « 500 nm). The com- 
puted flux distribution emerges around the minimum in the opacity 
caused by two collision induced absorptions, the H 2 -H 2 infrared vi- 
brational bands and the Lya wing broadening by H-H 2 encounters. 
Comparison in Fig. 1 10 | of our new syn thetic spectra (solid lines) 
with evaluations from Kowalski (2007) (symbols) for r c(f = 3000 
K and 5000 K shows an excellent agreement. 




Figure 11. [(B - V),(V - K)] and [(B - V),(V - /)] two-colour diagrams 
for hydrogen models with logg values of 7.0, 7.5, 8.0 and 8.5 (from top 
to bottom). The solid and dotted lines represent models with and without 
collision-induced wing of the Lya line, respectively. Select T e g values are 
labeled along the curves. Observations of DA (squares) and non-DA (trian- 
gles) WDs are taking from Bergeron et al. (2001). Evaluations from Kowal- 
ski & Saumon (2006) at log g = 8 and 7/ cff = 6000, 5000, 4000 and 3000 K 
are shown as open circles. 



The collision induced wing absorption does not affect signifi- 
cantly the structure of cool WD atmospheres since most of the flux 
is radiated at longer wavelengths. However, large deviations from 
the UV and blue colours are predicted with this opacity. In Fig.ll II 
we show the computed (B — V, V — K) and (B — V, V — I) two- 
colour diagrams for our new models (solid lines) compared to ob- 
servations (symbols) and previous model calculations without the 
Lya red wing (dotted lines). These atmosphere model sequences 
correspond to logg = 7, 7.5, 8 and 8.5 (from right to left). The 
effects of excluding the new opacity source are clearly visible on 
the B - V colour, giving rise to an increase of approximately 0.5 
mag. at the coolest models. This additional opacity does not affect 
the V and redder filters signifi cantly. The compar ison with the ob- 
served sequence of cool WDs ^Bergeron et al.l200lh shows that the 
new models reproduce the observations much better than the old 
models. In particular, we can see that the colours of the calculated 
models follow better the tendency of a linear sequence that the ob- 
servations indicate. On the other hand, the sequence at log g = 8 
displayed i n the (B - V, V - I) diagra m is very close to values 
obtained by Kowals ki & Saumo n (2006) at select and for the 
same surface gravity. 

Fig. [T2]displays a colour magnitude diagram using HST ACS 
filters (Vega-mag system) for DA- WD cooling sequences of sev- 
eral stellar masses in t he range . 52 < M(M ) < 0.88. A similar 
diagram was shown in iKowalskil feOOTl) for M = 0.5M G . Present 
calculations are based on a homogeneous set of evolutionary cool - 
ing tracks of hydrogen rich DA white dwarfs dRenedo et alj|2010h . 
These models take into account the most up-to-date physical inputs 
including the complete evolutionary history of progenitor stars, el- 
ement diffusion, chemical stratification, carbon-oxygen phase sep- 



Figure 12. (Mfg{4w,mf(,Q6 W - ffifgMty) colour-magnitude diagram for 
white dwarf sequences with masses 0.525, 0.570, 0.610, 0.660, 0.705, 
0.767, 0.837 and 0.878 M Q . Pure hydrogen atmospheres are assumed. Re- 
sults for selected T^g values are indicated with fill circles. Open squares 
correspond to the situation in which the Lya collision broadening is not 
considered for the 0.525 M Q model. 

aration, and crystallization processes. Fig. [12] shows that the Lya 
quasi-molecular opacity reduces the shift to the blue due to the H2 
collision-induced infrared absorption in the cool extreme of the se- 
quences. As a result of Ly-a opacity, the Mfgiw magnitude and the 
>nF606w - mF 8 i4iv color turn out to be 100-200 K cooler, thus imply- 
ing a larger age for a WD of a given stellar mass (see also Kowalski 
2007). 



5 CONCLUSIONS 

Kowalski & Saumo 3l l2006n (KS) identified the missing absorption 
in the blue and UV radiation of cool WDs as due to Lyman-a 
broadening by H-H 2 collisions. Here, we present independent cal- 
culations for the collision induced broadening of Lyman a line in 
dense gases typical of cool white dwarf atmospheres. Our results 
largely agree with those of KS. Quasi-molecular lines in these at- 
mospheres arise from radiative collisions of excited atomic hydro- 
gen with unexcited neutral molecules H2 or atoms H. As in the 
KS work, present line opacity evaluations are based on the quasi- 
static approach, which relies on the nearest perturber approxima- 
tion and the Franck-Condon principle. Line broadening processes 
are therefore described in terms of independent absorptions occur- 
ring during collisions between a hydrogen atom and a perturber, 
either another H atom or a H2 molecule. 

We revisit the properties of DA stars in the range 8000 K > 
r eff > 2500 K by analyzing the synthetic spectra with our im- 
proved models. The broadband colours located in the UV and blue 
spectral regions are shown to differ substantially from those pub- 
lished in previous studies which did not include collision induced 
Lya wing absorptions. Full tables containing both cross-sections of 
collisionally-broadened Lya wings and WD colours are available 
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at http://www.fcaglp.unlp.edu.ar/evolgroup or upon request to the 
authous at their e-mail addresses. 

The present study points out that a detailed knowledge of the 
simplest polyatomic molecule, H 3 , is of fundamental importance 
for a precise understanding of Lya opacity processes in cool white 
dwarfs. Accurate potential energy surfaces are required for their 
excited electronic states in order to precisely account for the effects 
of H-H2 collisions on the UV and visible spectrum of these stars. 
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